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We discuss the cosmological implications as well as possible observability of massive stable colored 
particles which often appear in the discussion of physics beyond the standard model. We show that if 
their masses are bigger than few hundred GeV and they saturate the halo densityand/or they occur 
, with closure density of the universe, they can be ruled out by existing limits on heavy stable particles 

■ from analysis of anomalously heavy isotopes of ordinary nuclei as well as from the observations in 
^\ ' WIMP detectors. We also derive constraints on their masses if they annihilate to produce gamma 

rays. We then comment on the possibility that these particles could be responsible for the ultra 
^jQ' high energy (UHE) cosmic rays with energies > 10^*^ eV observed in several recent experiments and 

I in particular point out that their low inelasticity argues against the possibility they can explain the 

■ observed UHE events. 
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Heavy stable particles carrying gauge charges naturally arise in many field theoretic contexts. Such particles with 
electric charges have been extensively discussed in the past ^ in connection with their possible role as dark matter. 
' We will assume in what follows that these particles are electrically neutral but carry color. Examples of such particles 
■ (to be denoted henceforth by X) can arise in supersymmetric models where the gluino ( rather than the photino ) 
\ is the lightest SUSY particle in which case we will have stable gg bound states with masses in the low TeV range 
(—! ■ (or much lower masses in the light gluino scenario @). There can be examples of gauge mediated supersymmetry 
Q_i breaking models where one can have color octet messengers that can be stable and heavy (with masses easily in 
' , the 10 TeV range). There are also examples of mirror models for particle where the particle spectrum as well as the 
gauge symmetries are completely doubled These models have been proposed in connection with understanding 
(— I \ the neutrino puzzles or to solve the strong CP problem. If in these models the color SU{'i)c groups of the mirror 
universe mixes with the that of the known universe, then the mirror quarks will share the strong interaction but 
not the observed electroweak interactions. Finally, the well known magnetic monopoles that arise in grand unified 
theories, although overall color neutral can carry a color cloud in addition to magnetic charge if the initial GUT group 
^ \ is simple. 

In this paper, we comment on the cosmological aspects and possible experimental observability of such particle^. 
We use a plausible ansatz for the annihilation as well as scattering cross-section for these particles and draw the 
following conclusions: (i) if their masses are above several hundred GeV, and they saturate the halo density, they 
should have been seen in the existing underground detectors and are excluded by present searches [Qjsj ; and secondly 
(ii) in the same mass range they are also excluded by the recent limits on anomalously heavy isotopes of ordinary 
nuclei j^]. Furthermore, if either the X particles are their own antiparticles (as in the case of color octets or gluino- 
gluon bound states) or the X and X are equally abundant ( saturating the halo density in both the cases), their 
annihilation in the halo can lead to energetic gamma ray fluxes. Present limits on these fluxes rule out the mass range 
Mx < 300 TeV. 

We also focus on recent suggestions in the literature to utilize particles [§,|oi in the above general category to 
explain the puzzling phenomena of ultra high energy (UHE) cosmic ray events p l[| observed in several independent 
experiments in the past decade. One could imagine these colored but electrically neutral particles |1Q] playing the 
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^For an earlier analysis of strongly interacting particles which constitute the halo dark matter see G. Starkman et al. j^. Our 
work is complementary to this and addresses several issues not covered in Ref. 



1 



role of the original high energy particle as the cause of these UHE cosmic rays. Our discussion leading to the above 
constraints also makes it clear that if these particles have masses above several hundred GeV, their hadronic collisions 
have very low inelasticity making them unsuitable as the UHE primaries which need to have significant inelasticities 
in order to generate the observed shower characteristics. 

This paper is organized as follows: in section 2, we present estimates of the annihilation, scattering and pair 
creation cross-sections of such particle. Using these estimates in sec. 3, we proceed to discuss the residual abundance 
of such particles in the present universe, which we find, could provide the closure density of the universe. In section 
4, we discuss their prospects for nuclear capture and formation of rare isotopes and their observability in various 
experiments. We then comment briefly on the connection between the UHE primaries and the massive X- particles 
in section 5. In section 6, we discuss in some detail the suggestion that the magnetic monopoles are responsible 
for the UHE cosmic ray events and argue that very low inelasticity of monopole-nuclear collisions tend to preclude 
this mechanism. 



II. REACTION CROSS-SECTIONS FOR THE X PARTICLES 



Because the X particles are much more massive than the usual hadrons, we must assume that their mass originates 
not from QCD but from physics at some new scale such as the scale responsible for breaking supersymmetry. Con- 
sequently we envision a small core for the X particle, of size 0{M^^) where the bulk of the mass is residing. This 
specific novel feature distinguishes the X particles from the ordinary hadrons. 

In X-nuclear or XX low momentum transfer collisions, only gluonic degrees of freedom residing in a hadronic cloud 
of size ^QC£, — Fermi need to be involved. Hence wc expect to have generic hadronic cross-sections 



crx-N 



KqId - 10 ™6 (1) 



The same estimate applies also to X — X elastic cross-section. On the other hand va X — X annihilation process, the 
cores of X and X particles should overlap and we therefore expect this cross-section to be much smaller: 

^^7^x « # « "'(10"' - 10"'°) G-eV-^ (2) 

X 

Likewise, just as for ti pair production, X — X production cross-section in hadronic collision is also very small: 

Mx 

^hadrons^XX ~ T7t/(^^) (3) 

with /(■^^) representing some further (structure function) suppression. 

The above crude estimates of the cross-sections for processes involving the X particles suffice for our discussion 
and unless the above intuitive arguments are completely off, the constraints obtained for the X particles will remain 
qualitatively unchanged. 

One may wonder whether the interactions of the soft hadronic clouds could modify the X — X annihilation cross- 
section. The soft hadronic interction extends over typical hadronic scale ( say, a Fermi) and also have typical hadronic 
strength i.e. effective potential depth of order C/q < GeV . If this interaction can pull together the "cores" of the 
two colliding X particles, then one would expect the cr^"Ijf to be of the same order as the Cx-x i-^- order 
o'hadron ~ {Fermi)'^ rather than the small value (Eq.2). The hadronic cloud interaction is strong enough to pull 
the "cores" of the X particles when the X particles are sufficiently slow moving, f3x < \J ^f^j with their center of 
mass kinetic energy less than a GeV. In this case we would expect a significant enhancement in their annihilation 
cross-section all the way to the value of say nucleon-antinuclcon annihilation at low velocities: 

^7-x « -^{Fermzf (4) 

This estimate will hold for example for relic X particles in the late cool stages of the universe. However the earlier 
estimate Eq. 2 will hold in the early hot stage when the center of mass kinetic energy of the collision exceeds 
Uq « ^QCD and we can ignore the cloud effect. We will see in the next section that this result has a profound 
implication for the number density of the relic X particles and their subsequent behaviour. 

As we will see in the next section, the freezeout temperature T* depends only logarithmically on 0"^'",^ . The ratio 
T* /Mx is always larger than 0.02 so long as cr""" < (Fermi)'^. For the value of Mx ~ 1 — 10 TeV (as is the case we 
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are considering), the temperature at decoupling T* « 20 — 200 GeV, is considerably larger than Aqcd making the 
enhancement of the cross-section noted above irrelevant for the relic abundance. 



III. RELIC DENSITY OF X PARTICLES 



In the very universe T 3> Mx , we expect the density of X particles to be comparable to that of other relativistic 
particles. However as the universe cools below Mx, the annihilation of the X — X pairs dominates over their 
production gradually decreasing the ratio ux/tLj, until the freezeout temperature where rate of annihilation becomes 

slower than the expansion rate of the universe. The annihilation rate is given by ^^^^ = a^^vnx so that the 
freezeout temperature is dictated by the following inequality JTst : 

«s < i.66V^r^ 

ttM^ - Mpi ^ ' 

Using next nx ~ gx( ^2"^"^ )'^^^^~^^ ' (where gx is the number of degrees of freedom of the X particle: e.g. for a color 
octet spin zero field, gx = 8 etc) we can rewrite (Eq. ) in terms of the dimensionless quantity ^ = as follows: 

(2^)5/2 W 



Taking Mx/Mpi = IQ-^^ - 10-l^ g* ~ 10 and ~ .12, we find 

'(27r)5/ 

Finally using « g*T^ we have 



e « 1/2^ + /nlO^^ + « 30 (7) 



- '-''^^ « 170e^ « 5 X (10-- - 10--) (8) 



2A^arMpi ^ Mpi 

Using next rjB = — ~ 6 x 10"^°, we find 



Hence the mass densities are in the ratio 

px_ _ nxMx 
Pb ripmB 



^ ~ 10-2 - 10-3 (9) 
ub 



= 1-100 iforMx -1-10 TeV) (10) 



Thus we find the nice feature that taking AIx in the range of 1-10 TeV suggested by independent considerations 
naturally results in closure density of the universe. This is of course not an entirely new result but clearly related 
to the discussion of the lightest supersymmetric particle (LSP) with a corresponding weak annihilation cross-section 
being the dark matter of the universe. 

We will now show that there are severe difficulties with such relic particles if the X have hadronic strength interac- 
tions with ordinary matter. Once the temperature T < Aqcd — -3 GeV , the QCD phase transition temperature, the 
gluon and quark gas disappear. At this hadronic recombination temperature, the purely strogly interacting particles 
almost decouple from the baryons (and from the rest of the background radiation). Once px > P7, a standard cold 
dark matter scenario may be implemented: the X particle would start forming structure and galactic halos in par- 
ticular. Subsequently the potential wells of these particles will trap the dissipating baryons which form the galactic 
discs. If we assume that the local halo density p w 0.3 GeV/(cm)'^, is generated by our X particles, we will have 

nx =i 0.3 GeV/Mx =i 10-^ cm'^ (11) 

This of course significantly higher than the value of relic density mentioned earlier (using w 300 cm~^ one finds 
that n^^'*'^ ~ 10-^" — 10""'^^). We will now discuss three possible manifestations of the X- particles: (i) in underground 
detectors ; (ii) via rare isotope formation and (iii) via annihilation in the halo. This discussion will be done for both 
cases of galactically clustered as well as cosmological, unclustercd, X particles. 
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IV. MANIFESTATION OF THE X PARTICLES 



A. Wimp detectors and X particles 



Let us first consider issue (i) - namely the possible manifestation of the X particles in underground detectors if the 
X's saturate the halo density. Even though the X particles are strongly interacting (unlike the usual wimp candidates 
for dark matter) their penetration depth in earth or water can vastly exceed that of the normal hadrons. This is 
so despite the fact that as indicated in sec. II, the typical X-nucleus scattering may well be of the size of normal 
hadronic cross-section- i.e. a « ttR^^A^Z) with R ~ 1.2A^/^ Fermi, the nuclear radius. The reason for this is the 
large mass mismatch between the X particle and a normal hadron. This implies that we require many collisions 

Mx 

ncoii ~ Mx/ < m(A,z) >- (12) 

(where we have chosen < m(^_2) >« 20 GeV for average water-air-crust nuclei) before the initial X particle appre- 
ciably slows down or is deflected from its path. The effective penetration depth is then given by 

tp = Xricou (13) 

where A is the mean free path between succesive collisions given by 

1 ^ lO'^ cm 

~ n{A, Z)a{X - {A, Z)) ^ a/{FermiY ^ ' 

where we used n{A, Z) = NAvagadroP/ < m{A, Z) > and p sa 3 gm cm~^, the average density of the crust of the 
earth. As a result the penetration length is given by £p = A 2o*cfeV ~ ^/(FermV ( 20^^GeV ) • Thus for Mx > 2 TeV, 
the penetration depth is about a kilometer which is the typical depth for most underground Wimp detectors. Let us 
emphasize again that the key to the large number of the collisions, leading to large penetration depth is the small 
energy loss in any given X-nucleus collision due to the large mass of the X particle. Even if the collision "drags" 
the nucleus to the same velocity Px of the X particle, the kinetic energy of the nucleus will only be ^m{A, Z)0^ 
i.e. a factor m{A,Z)/Mx smaller than the kinetic energy of the X particle. Hence a large number of collisions 
ricou = Mx/m{A, Z) will be required to reduce its energy by a factor e (i.e. (1 — j^^)^""" = e~"^). 

A somewhat more careful estimate of the X— nucleus cross-section cr{X — {A, Z)) can be obtained as follows. It has 
been argued before |lj] that the slow XG (where G stands for the gluon) color neutral bound state, to which we will 
continue to refer to as the X particle, sees the nucleus as a potential well of radius R{A) and depth of order Vb — 10 
MeV. The scattering is predominantly in S-wave. For the relative motion with reduced mass m ~ m{A^ Z), energy 

« "^^^'^^f^ ~ 10 keV and momentum p ~ to(A, Z)/3 w 20 MeV (assuming virial velocities for the particles) one gets 
a phase shift 



-pR + arctan tan\/ p^ + 2mVoR I w —pR (15) 



where we have used the fact that \j2jmVQR 3> 1 so that the cross- section is essentially geometric i.e. 

a ^t:R^ 7=iTxA^I^{Fermif (16) 

Using this one gets via Eq. 13, 14 and 16 a penetration depth « 100 meters in the rock for Mx ~ 2 TeV, confirming 
the earlier estimate. 

If the above analysis applies, then the "wimp" detectors would be able to detect the full flux of the X particles if 
located at appropriate depths where the enrgy attenuation is very small. Whereas most underground detectors (e.g. 
lH) are located at greater depths, some are relatively shallow ( 0]) and will suited for this search. Note further that 
since the X— nuclear cross-sections are about k, 10^^ times bigger than the generic heavy Dirac neutrino cross-section, 
if the local X particle density is either given by the estimated halo density above or even the much smaller relic 
density, the detectors that ruled out the heavy Dirac neutrinos would a-fortiori also rule out the heavy X particle. 

We should note that the X particles can be captured in nuclei once the center of mass collision kinetic energy 

\Mx0^ is smaller than the expected few MeV X-nuclear binding energy i.e. /3x < \ ^° ~ which is satisfied 



i^ivixH niiictiici tiidii Liic CAjjCL.tcu icw ivic v yv-iiuL.iccii uiiiuiiig ciitsigj i.e. MA y m{A,Z) 50 

for virial velocities. However once the nucleus is captured onto the X particle, its crosssection (hence its attenuation 
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length) will not necessarily increase (decrease) and in principle the nucleus could be dislodged in subsequent collisions. 
Hence any possible capture of the X particle need not change the above estimate. However accompanying the capture 
of nucleus by the X, there could be some rather striking features as 7 or neutron emission which could serve as more 
sensitive signatures although it must be remembered that extra photon emission rate will down by a factor of aem- 



B. Rare X-isotope formation 

We have estimated above that the relic cosmological abundance of X-paritcles relative to hydrogen is about 10""^ 
or so. Clearly this ratio need not be maintained locally on earth where we have a much enhanced baryon number 
density. To estimate the terrestrial X particle effect, we first note that a flux of such particles impinging the 
earth thru Hubble time yields a total of ^x-tn — 2 x 10^^ ^x X -particles falling on each square centimeter on the 
earth. These X particle would most likely accumulate in the ocean. For an ocean depth ^ 10 kilometers, and surface 
hydrogen density of nn — jgNAvagadroW^ = lO'^^H/cm'^, th = ^ ^ 10~^^^x and a similar number obtains for 
oxygen also. After slowing down to thermal velocities the X-particle will be captured by the nuclei in the ocean 
water- most likely the oxygen nuclei- forming anomalous heavy (« Mx) isotopes. A careful search for precisely such 
heavy isotopes have been carried out in a recent experiment by Hemmick et al |q| and they find that the upper bound 
on the ratio roxygen < 4 x 10"" - 3 x 10"" and th < 2 x 10'^^ - 3 x 10^'^ (for Mx = 0.1 - 10 TeV). Even if 
the collection period is shortened to 100 million years (a reasonable lifetime of the ocean), these upper bounds would 
seem to exclude X particle fluxes $x > 1 cm^^ sec^^. If the AT's constitute the dark galactic halo the expected flux 
$jsf can be estimated to be 

$x « nxvx = 3 X 10^ ~ 10^ - lO'' for Mx = 10 - 1 TeV (17) 

Mx/ GeV J ^ \ J 

Thus AT as a halo dark matter is ruled out by the experiment in Ref. ||^. On the other hand, if AT-particles do not 
dominate Pcosmos and do not cluster in the Halos, so that Eq. 8 apphes, then $x « nxVx ~ 6 x (10^^ — 10^^) which 
given the crudeness of our estimate is marginally consistent with the above lower bound. 



C. Annihilation in tlie Halo 



In this section, we consider the possibility that if the X particle is its own antiparticle or the X and X are present 
in the halo in equal abundance, their annihilation can give rise to energetic 7-ray fluxes. Present bounds on such 
fluxes can then lead to constraints on the X particle density nx- We derive this below. 

We pointed out in section 2 that at small velocities f3 « 10"'^, where the X energies are below Aqcd, the X — X 
annihilation is enhanced and is given by axx ~ -^{FermiY . The X particles in the halo will then annihilate at a 
rate 

= nxvxc^xx ~ nx(Ferrm) c (18) 

nx 

with nx ~ Mx /GeV ^ 10"'' cm"^ for Mx = 3 TeV. This implies a halo lifetime of 10^° sec, considerably exceeding 
tn- As a result, only a fraction « 10"'^ of the initial X's will annihilate in Hubble time. Thus any galaxy will lose 
w 10^'^ of its rest mass via such annihilation. This exceeds the total energy output of stars. 

The problem however is much more severe. The stellar radiation is largely carried away by the optical photons, 
whereas roughly 40 per cent of the energy released in X — X annihilation will convert to 100 GeV to TeV photons. 
There are extremely strong bounds on these fluxes on earth [ p^ : 

$^ < 10"^ cm~'^ sec^^ sr^^ (19) 

Since ~ nxcrxx'^^xL, Eq. 9 implies nx < 10^® cm~^ where we have used L ~ 9 kilopersecs as a typical halo 
radius. 

In addition, these decays will also inject a large number of relativistic e+e" which if trapped in the galaxy will have an 
energy density of ~ 0.3/ GeV cm~'^ where / ~ lO^'^ is the decaying fraction. Using equipartition of energy, this should 
not exceed the energy density in the galactic field which can be estimated to be ub — ^^^^^ ^Gauss) q 2 ^.j^^-s^ 
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However, a naive estimate of Ug+e- gives « 3 x 10~* GeV ~ 3 x 10^ eV which is 10^ times larger. Thus both these 
arguments would imply that the halo density cannot be dominated by the massive X particles. For unclustered 
X-particles, the analog of Eq. 19 with L the effective halo size replaced by the Hubble radius, we find nx < 10~^ 
cm"'^ (ignoring small redshift effects). Note that this is comparable to the cosmological prediction for nx in Eq. 8. 



V. MASSIVE COLORED PARTICLES AS THE SOURCE OF UHE COSMIC RAYS 

During the past decades cosmic ray events with energies greater than 10^° eV have been observed in several 
independent experiments. This may call for new physics since cosmic ray proton with E > 5 x 10^ eV, scatter 
off the cosmic microwave background producing nucleon resonances leading to the Griesen-Zatsepin-Kuzmin cutoff. 
Moreover, the mean free path for this interaction is about 6 Mpc (or redshift z — 0.0f25). Since such energetic protons 
are unlikely to be bent by the galactic magnetic fields their direction will point directly to the source which must be 
less than 6 Mpc in distance if protons are to be the source of these UHE showers. No such source has been found 
I p^ . On the other hand in the same general direction, there are some viable sources at greater distances (e.g. the 
quasar 3C 147 at a distance of 240 Mpc or the Seyfert galaxy MCG 8-f 1-f 1 at a distance of about 10 Mpc). For the 
primaries to be coming from these distant sources, they must avoid the analog of the GZK cutoff and must be strongly 
interacting to produce showers [|o|. Therefore electrically neutral, colored particles are natural candidates for the 
primary particles. Furthermore, in order to successfully reproduce the observed extensive showers, these particles 
should have inelasticities comparable to those of the usual primary cosmic rays, namely protons or nuclei. In primary 
proton collisions, about half of the initial energy is utilized to make the energetic secondaries. The latter inelasticity 
is high. The arguments of the previous section implies that if the UHE primaries were heavy (> few hundred GeV) 
neutral colored particles they would have been seen in the underground detectors. Of course, if they are light (as in 
the model of Chung et al. Q), our arguments do not apply. We now turn to the suggestion that heavy monopolcs 
of grandunified theories may explain the UHE cosmic ray showers 113]. 



VI. MONOPOLES AS SOURCE OF ULTRA HIGH ENERGY COSMIC RAYS 

It was suggested in Ref. that monopoles of grand unified theories may be responsible for the ultra high energy 
(UHE) cosmic rays with energies above 10^" eV II4]. The attractiveness of the suggestion arises from the fact that 
this offers a natural mechanism for acceleration of the particles to such ultra high energies by the galactic magnetic 
fields. For natural values of the monopole charges expected in GUT theories i.e. Em ~ VNgMBcLcoh — 10^° eV 
where = qe/'ia is the monopole magnetic charge, N is the number of coherent domains traversed and Lcoh is the 
typical length of a coherence domain expected to be of order 100 persecs. To mimic the observed cosmic ray events, 
the monopole mass tum must be less than 10^° GeV since it must be relativistic. There is also a kinematically a 
maximal value for the inelasticity i.e. energy transferred from the monopole in the scattering process on a target. To 
obtain this the authors of Ref. 7 note that the maximum energy transferred in a collision of a monopole of mass niM 
and energy Em on a target nucleon of mass m is given by: 

2mEM 

Qmax — JZ 7; ' i 

[2mEM + mfi) 

obtained for backward center of mass elastic scattering. We see that for UHE case of Em — 10^^ GeV, C acquires its 
maximal value (C = 1/2) for tom '^-^ 10^ GeV and decreases for higher values quadratically. If this maximum value 
for the inelasticity were attained, it would be very favorable for the interpretation of the monopoles as the source of 
UHE's. To see if this maximal inelasticity value can in fact obtain, one has to consider the dynamics of the monopole 
collision. We find that such high values for the inelasticity are unlikely to be realized in practice if either one of 
the two collision mechanisms that we consider dominates. Below we describe these two monopole-nucleus collision 
mechanisms and estimate the reaction rates and the inelasticities in both cases. 



A. Hard electromagnetic scattering of monopoles 



The most obvious reaction mechanism for initiating the UHE events by the monopoles utilizes their magnetic 



interaction with charged particles. The dominant mechanism briefly alluded to in Ref. [ 
due to the Lorentz boost generated large transverse i?-field seen by the stationary nuc 



12| involve nuclear dusruption 
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with 5 being the impact parameter and 7 — = lO^'^^/mM iJT. GeV . This field can disrupt the nucleus by accelerating 
the protons (relative to the neutrons). Let us estimate how big this effect is. 
The external perturbation seen by the stationary nucleus is: 

H'{t) = cEt ■ Sri/(i) (22) 

where ri is the ith proton's location relative to the center of the nucleus and f{t) is a pulse shape normalized as 

f(t]dt ~ At ~ — (23) 
7c 

where b is the impact parameter and the 1 /7 factor is due to squeezing towards the transverse direction of the B and 
E field lines of the moving monopole. This leads to a 7 independent product EAt. This product in turn controls both 
the boost invariant transverse momentum imparted ~ FAt and the excitation probability. The probability of exciting 
the nucleus is given in first order time dependent perturbation theory by the following sum over all non-ground states: 

P = SPo^n < egM^'S'S„^i J |/(t)e'^^"*dtn < 0|ri|n > ^ (24) 
< (ZEfl J f{t)dt\^ < 0\r^\0 >w Z^iEAtfR" (25) 

« ^2i?LcJ6' (26) 



where Schwarz inequalities, completeness and eqM = 1 have been used. Thus only if 6 < -^Y^nuci ~ 

Z X 1.2A-^/^ Fermi « 60 Fermi, will we have P ~ 1. The excitation cross-section is therefore smaller than 
CTdis ~ TT^^ ~ cxr? . Furthermore, in each disruption, the total energy loss in the excited system rest frame 

does not exceed the total nuclear binding energy i.e. AE* < 100 MeV. This gets boosted at most to AE = jAE* in 
the laboratory frame. Even in this extreme case the fraction of the initial monopole energy lost in GeV in a collision 
is 

V < {jmM)~\jAE*) = AE*/mM - 0.1 GeV/mM (27) 

Along the line of sight in the atmosphere of average length £ ^ 20 km, there are at most Nc — £ rLMitrogen<^dis < 10^ 
such collisions and the total energy loss is therefore Ncrj ~ lO^/mA/, which is very small indeed for values of the 
monopole mass of interest. 

Let us next consider the "hard collision" of a point like monopole with a point like charge. For a classical collision at 
impact parameter 6, the boost invariant transverse momentum transfer is given by ApT — eExAt « ^p^6 — equ/b, 
where as discussed above Et = Bt^ — equ^ jb^ and At k, b/^. This is precisely the same as the momentum transfer 
in a hypothetical Rutherford scattering with charges ^ Zi = e and Z2 = qM- Carrying this analogy to the field 
theoretic domain, we write down the Rutherford like cross-section for magnetic monopole and charge scattering as: 



dt ''"'"^^ - 



,„,e(0 ^ (28) 



Although formally this leads to a divergent total cross-section, the invariant transfer in our case is simply 2mT where 
m is the mass of the electrically charged particle and T is the laboratory energy transferred to it. The effective 
cross-section for monopole scattering obtained by weighing ^ with the inelasticity T/Em — t/2mEM is given by 



^ It is amusing to recall that in nonrelativistic classical mechanics there is an exact description of the motion of a point charge 
in the field of a monopole. It is a convergent and a reflected spiral restricted to the Poincare cone with vertex at the origin 
and opening angle tt — Oq with Oq eg / L = eg/pb, the ratio of the electromagnetic to the orbital angular momentum, such 
that when the cone is opened into a plane, the spiral becomes a zigzaging straight line The actual final scattering angle 
is 6 ~ 60 and is just the same as for a Coulomb problem with Q = gM at the origin. Independently of all the above, we can 
compare the scattering amplitudes for charge-charge (e — e), monopole- monopole [qm — qm) and charge -monopole scattering 
(e - qm) via the inequality Ae-e Aq^jq^, > Al_q^j. 
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^^^^ = ^ / i^T^'' " ^^"(^™-/^-) (29) 

For i?M of order 10^^ GeV, this gives a cross-section for monopole scattering of magnitude 10^**° cm^, which is utterly 
neghgible. 



B. Nuclear interaction of monopoles 



The magnetic monopole could have a hadronic cloud consisting of gluons and qq pairs. This can happen for example 
if the gauge group whose breaking leads to the monopoles includes SU (3)c as a subgroup. The hadronic cloud will 
then have a typical dimension ~ ^qcd order one Fermi or so. The hadronic cloud is inherent to the monopole and 
is expected to regenerate after every collision. Thus the monopole will have a hadronic cross-section for interaction 
with nuclei similar to that of a proton. This feature would seem to favor this scenario for understang the UHE events. 
In particular, one event from the Yakutsk experiment appears to be neutron rich suggesting a hadronic origin for the 
UHE's. 

Unfortunately there is a crucial difference between the inelasticities of a monopole initiated and proton initiated UHE 
reaction. The hadronic outer layer of the monopole carries only a minisculc fraction e = GeV/niM of the monopole 
mass. This can be seen for the 't Hooft-Polyakov monopole, where the bulk of the monopole mass Mm ~ -^g/sg 
(where Mq and go are the symmetry breaking scales and the gauge couplings associated with the GUT group whose 
breaking causes the appearance of the monopoles) resides in the innermost regions of the monopole (of dimension 
Mq^) and an order TeV contribution from the weak layer of the scale dimension . The extended hadronic layer 
contributes only a fraction Kqcd I Qs ~ GeV to the monopole mass. In the language of the parton model, the gluons 
and the qq pairs are the "wee" partons in the very small x- region of the parton distribution in the monopole. Since 
only thses hadronic constituents can participate in the hadronic monopole nucleus collision, we expect the inelasticity 
to be of order C < 2; < GeV/niM- This result is clearly in contrast with the characteristic inelasticities of order 1/2 
in proton nucleus collision. 

Possible binding of the monopole to the nuclei is also not likely to be important. The reason for this is that the 
binding potential is given by 

^^-^9M (30) 

where fj, is the magnetic moment of the nucleus and is the magnetic charge of the monopole. Due to relatively small 
magnetic moment of the Oxygen and Nitrogen nuclei, it is not clear whether these nuclei can sustain the requisite 
p-wave monopole-nucleus bound state |Q. Even if such a bound state formed, it would strip off in the first nuclear 
collision and the rate of monopole nuclear capture to reform the composite is likely to be slower than that of the 
nuclear reactions destroying it. Hence the mechanism of enhancing nuclear monopole interaction by binding them to 
nuclei and thus enhancing xh from GeV/mm to xh ~ A(GeV)/mM would be rather inefective. It turns out also 
that endowing the monopoles with large cross-sections to catalyze baryonic decays a-la Callan and Rubakov does not 
help either^ 

Thus in summary, we have shown that massive, stable and neutral colored particles with masses more than a few 
hundred GeV are ruled out by exisiting WIMP searches and masses in the range 100 GeV to few TeV by the present 
results on searches for anomalous heavy isotopes. They are also not suitable for explanation of the UHE cosmic rays. 
Furthermore, as far as the monopoles are concerned, while we have not conclusively established that their interactions 
cannot lead to the UHE's, we have raised some plausible doubts about the effectiveness of the mechanism due to very 
small inelasticities in monopole nuclear interaction which can be inferred from simple intuitive arguments. To clearly 



^ While this general issue has been discussed in the past (see e.g. |19|), the following argument may still be worth presenting. 
Let us consider a blue or red giant star of radius i? « 3 x (10^^ — 10^'') cm. and column density of « ~ 10* — 10^ gr / cm^ . 

Clearly if the monopoles lose even 10""^ of their energy in our atmosphere, they would stop upon hitting these objects. During 
its lifetime of 10 million years, the number of monopoles accumulating there is Nm ~ Att^^mR^T « $m10''^^^. After the 
gravitational collapse to a neutron star, most of these monopoles will reside in the nuclear core catalyzing one nuclar decay in 
each hadronic time scale of 10~^^ sec. Altogether in 3 x 10^® sec, each monopole catalyzes 3 x 10^® decays destroying a 
fraction 10~^® of the core. To allow the neutron stars to live for a billion years, we must require that Nm < 10^* implying that 
$A/ < 10^^'**^ cm~^ sec~^ sr~^ . This would totally exclude the above scenario. 
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confirm or rule out the mechanism, detailed simulation of the effect of relativistic monopole air nuclei collision in the 
various experimental setups is required. Our net conclusion is that the monopole scenario for the UHE's is rather 
unlikely. 
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